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Abstract 

We study the modification by the presence of a plane wall of energy level shifts of two-level 
atoms which are in multipolar coupling with quantized electromagnetic fields in a thermal bath in 
a formalism which separates the contributions of thermal fluctuations and radiation reaction and 
allows a distinct treatment to atoms in the ground and excited states. The position dependent 
energy shifts give rise to an induced force acting on the atoms. We are able to identify three 
different regimes where the force shows distinct features and examine, in all regimes, the behaviors 
of this force in both the low temperature limit and the high temperature limit for both the ground 
state and excited state atoms, thus providing some physical insights into the atom-wall interaction 
at finite temperature. In particular, we show that both the magnitude and the direction of the force 
acting on an atom may have a clear dependence on atomic the polarization directions. In certain 
cases, a change of relative ratio of polarizations in different directions may result in a change of 
direction of the force. 
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I. INTRODUCTION 



Quantum fields constrained by the presence of boundaries exliibit many interesting prop- 
erties because tlie presence of boundaries disturbs tlie modes of quantum fields. As a re- 
sult, fluctuations of quantum fields are altered, giving rise to a lot of novel effects, such 
as the Casimir [ll and Casimir- Polder forces j^, the light-cone fluctuations when gravity is 
;he Brownian (random) motion of test particles in an electromagnetic 



quantized 
vacuum 




and modifications of radiative properties of atoms in cavities 



such as the natural lifetimes and energy level shifts which have begun to be demonstrated 
in experiment [121, 13] • 

In the present paper, we are particularly interested in the Casmir-Polder force at finite 
temperature. The Casimir force, which has now been measured with precision [l^, is 
not only fascinating in theoretical research but also increasingly important in technological 
applications, such as nanophysics [l6|, chemical identification of surface atoms via atomic 



force microscope 



171] and the construction of novel biomimetic dry adhesives 18|. 



Casimir-Polder force at zero temperature is usually regarded as a result of the reshaping of 
vacuum fluctuations induced by the presence of boundaries. In a thorough understanding of 
the Casimir-Polder force one must also consider the effects of thermal fluctuations in addition 
to the vacuum ones. It should be pointed out that in the studies of thermal Casimir-Polder 
brce, a majority of researches invokes a macroscopic approach to atoms pioneered by Lifshitz 



19 



20 



2l| or a linear response theory of atoms 22|, |23| , although some uses a QED treatment 



of the atom-field coupling 



24 



25 



26|. 



With an atom-field coupling, the Casimir-Polder force could be understood as a result of 
the position dependent modification of radiative properties of atoms (the energy level shifts 
in particular) in the confined space in QED quantitatively. Let us note, however, that a clear 
and simple physical interpretation does not always emerge from the QED description. One 
interpretation is that the boundary induced effects on atom's radiative properties may be 



exp 



ained by the modification of vacuum fluctuations of the quantum electromagnetic field 



28 



29 



30| . On the other hand, the same effect can also be described as a consequence of 



2 



the reaction of the instantaneous atomic dipole to its own radiation field refiected from the 



boundary 



31 



,132, 



33 



3J]. In fact, it has been shown within QED that the extent to which 



each mechanism contributes to the total effect can be chosen arbitrarily just 



35 



by; 



36 



changing; 



37|, 



38|. 



the ordering of the atom and field operators in the interaction Hamiltonian 

As a result, there exists an indetermination in the separation of effects of vacuum ffuc- 
tuations and radiation reaction such that distinct contributions of vacuum ffuctuations and 
radiation reaction to the radiative properties of atoms do not possess an independent phys- 
ical meaning. Therefore, although quantitative results for spontaneous emission and energy 
level shifts are well-established, the physical interpretations remained controversial until Dal- 
ibard, Dupont-Roc and Cohen- Tannoudji(DDC) argued in 39| and 40| that there exists a 
symmetric operator ordering of atom and field variables where the distinct contributions of 
vacuum fiuctuations and radiation reaction to the rate of change of an atomic observable 



are separately Hermitian. Recently, this 
tive properties of atoms in various cases 



brmalism has been employed to study the radia- 
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42 



43 



44 



45 



46 



47 



48 



49 



51 



52|. In 



the present paper, we would like to use the DDC formalism to calculate the modification 
of energy shifts of atoms by the presence of a boundary in a thermal bath and the forces 
thus induced between the atom and the boundary. This is not a new problem. However, by 
using the DDC formalism which separates the contributions of thermal fiuctuations (includ- 
ing vacuum fiuctuations) and radiation reaction and allows a distinct microscopic treatment 
to atoms in the ground and excited states, in contrast to the macroscopic approach where 
atoms are treated as a limiting case of a dielectric, we hope to gain some new physical insight 
into the atom-wall interactions. Let us note that the same problem at zero temperature has 
already been discussed in using the DDC formalism 4l|. 

The paper is or gan ized as follows, we give, in Sec. II, a review of the general formalism 



developed in Ref. 



401 1 and generalized in Ref. 43|, |45| to the case of a neutral polarizable 



two-level atom interacting with quantized electromagnetic fields at finite temperature T, 
and apply it, in Sec. Ill, to calculate the position dependent energy shifts, which give rise 
to an induced force acting on the atom. Then we will study, in detail, this force in both the 
low temperature limit and the high temperature limit for both the ground state and excited 



3 



state atoms. Finally, we will conclude in Sec. IV with a summary of results obtained. 



II. THE THERMAL FLUCTUATIONS AND RADIATION REACTION CONTRI- 
BUTIONS 

We consider a neutral polarizable two-level atom in interaction with quantized electro- 
magnetic fields in a thermal bath of temperature T in the presence of an infinite conducting 
plane wall located at 2; = 0. We assume that the atom is at rest at a distance, from the 
wall, and has stationary states |— ) and |+), with energies —\tvjjQ and +\tujjQ and a level 
spacing tvjjQ. 

To investigate the modification of the atom's energy shifts by the presence of the plane 
boundary, we will first review the general formalism developed by DDC which allows a 
distinct separation of the contributions of thermal fiuctuations (including zero temperature 
vacuum fiuctuations) and radiation reaction to the energy shifts of atomic levels. 

The atom is supposed to interact with the quantized field in the multipolar coupling 



scheme 531], so that the Hamiltonian that describes the atom-field system with respect to 
the proper time r can be written as 

H{r) = HA{r) + Hp{t) + Hi{t) , (1) 

where 

HA{r) = hwoS^ir) , (2) 

k 

Hjir) = -^(r) ■ E(a;(r)) = • E(x(r))[^+(r) + S.ir)] • (4) 

Here k denotes the wave vector and polarization of the field modes, /x the atomic electric 
dipole moment, and the electric field operator, E, is evaluated along the trajectory x{t) 
of the atom. Sz, 5*+, S- are the pseudospin operators of the two-level atom, and ^^(O) = 
- S^O) = |+)(-|, S_(0) = |-)(+|. 
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With the Hamiltonian of the system given, we can write down the Heisenberg equations 
of motion for the dynamical variables of the atom and field from Eq. ([1]). The solutions of 
the equations of motion can be split into the two parts: a free part, which is present even 
in the absence of the coupling, and a source part, which is induced by the interaction of the 
atom and the field. We assume that the electromagnetic field is in a thermal state with the 
density matrix p = Q-PHp/{nc) ^ where /3 = he/ {TkB) is the thermal wavelength, and the atom 
is initially in the state \b). ks is the Boltzmann constant. To identify the contributions of 
thermal fluctuations and radiation reaction to the radiative atomic energy shifts, we choose, 
following DDC, a symmetric ordering between the atom and the field variables and consider 
the effects of E-l" (corresponding to the effect of thermal fluctuations) and E'"^ (corresponding 
to the effect of radiation reaction) separately in the Heisenberg equations of an arbit rary 
atomic observable G. Following the procedures that have been shown in Refs. 0, Q, 45], 
we take the thermal average over the electromagnetic field and decompose various rates of G 
into a Hamiltonian part, describing the energy levels of the atom shifted due to coupling with 
the electromagnetic field, and a non-Hamiltonian part, describing the effects of relaxation, 
then in a perturbation treatment up to order /x^, we find 
'dG{ 



— ) P) = Tl^tfrri^)^ ^i^)] + non-Hamiltonian terms , (5) 

/tf,rr / n J' 



where 



<^(r) = / rfr'(C5),(x(r),x(r'))[^f(r),/if(r')] , (6) 

H^i^ir) = -Yhj dr\x!,)M'^U{r')){^^i{r\^^^,{r')] . (7) 

Here [ , ] and { , } denote the commutator and anticommutator respectively while subscript " 
tf " stands for thermal fluctuations and "rr" for radiation reaction. The statistical functions 
{C^)i3 and {xfj)i3 are defined as 

(C5),(x(r),x(r')) = liiElixir)), Elixir'))}), , (8) 
(x5)/3(x(r),a:(r')) = ^([i?/(x(r)), E/(x(r'))]), , (9) 



which are also called the symmetric correlation function and the linear susceptibility of the 
field. Taking the expectation value of Eqs. and ([7]) in the atom's initial state \b), we can 
obtain the contributions of thermal fluctuations and radiation reaction to the energy shifts 
of the atom's level \b), 



i^E,h = ~l rfr'(Cj)^(x(r),x(r'))(x3Mr,r'), (10) 



{5E,\r = -\l dT'{xfi),{x{r),x{T')){Cf,U^,T') , (11) 



TO 



where {C^)b and {xfj)b, the symmetric correlation function and the linear susceptibility of 
the atom, are defined as 

(C^^Mr,r') = i(fe|{/if(r),/.J(rO}|6), (12) 

(X^).(r,r0 = ^(6|[/xf(r),/zj(r')]|6), (13) 

which are characterized by the atom itself. Explicitly, the statistical functions of the atom 
can be given as 

iCf,W,r') = i$^[(fe|/i.(0)|d)(rf|^,(0)|6)e^--(— ') + (6|/.,(0)|rf)(rf|/^.(0)|6)e--(— ')] , 

d 

(14) 

ixiW,r') = ^EK&l/^^(0)M)(rf|/^.(0)|6)e^--(— ') - (6|/i,(0)M)(d|/i.(0)|6)e--(— ')] , 

d 



(15) 

where um = Uf, — Ud and the sum extends over a complete set of atomic states. In order to 
calculate the statistical functions of the electric field, Eqs. ([H]) and we will use the two 
point function of the four potential, A^{x), at finite temperature, which can be obtained by 
the method of images both in imaginary-time and in space. In the presence of a boundary, 
one finds in the Feynman gauge 9j, 

D^.'^ix, x') = {A^{x)A'^{x')), = D;^}_(x - x') + D^;,^,{x, x') , (16) 



where 



PJreeKX ^ ) 4^2^^^ ^Z^^ [^^^ _ ^ -^^ _ -^^2 _ (3; _ ^/)2 _ _ ^/)2 _ _ ^/)2] ' 

(17) 

is the two point function in the free space and 



^TT^^oc \{tc - t'c + - isY - {x - x'Y - {y - y'Y - {z + z'Y\ 

(18) 

represents the correction induced by the presence of the boundary. Here Eq is the vacuum 
dielectric constant, the subscript " hnd " stands for the part induced by the presence of a 
boundary, e +0, ?7'^'^=diag(l, — 1, — 1, — 1) and the unit normal vector = (0,0,0,1). 
Then, we can write the electric field two point function as 

{EMr))EMr')))p = {EMr))EMr')))^Jree + {EMr))E,{x{r')))pMd , (19) 

where 

{EMr))EMr')))(Sjree = ^ E 

k=—oo 

1 

^ (x - x'y + {y- y'f + {z- z'f - {tc - t'c + ikl3 - ief ' 

(20) 



and 



a_ OO 

he 



{Ei{x{T))Ej{x{T'))),3^bnd = E [{^ij -'^ninj)dod'Q- did'j' 

k=—oo 



1 

X ■ 



{x - x'y + {y- y'f + {z + z'f - {tc - t'c + ikl3 - ief 

(21) 

Here &j denotes the differentiation with respect to x'y So, the statistical functions of the 
electric field can be obtained from Eq. (fT9l) as a sum of the free space part and the boundary- 
dependent part. 



Our aim in this paper is to study the modification of atomic energy shifts by the presence 
of the conducting plane boundary in a thermal bath and the force thus induced on the atom. 
So in Eqs. fllOl) and fill I) , the contributions of thermal fluctuations and radiation reaction to 
the energy shifts of the atom's level, we are only interested in the boundary-dependent part, 

iSE,)1f = ~jyr'{Cf^),,U^{T),x{T')){xiUT,T') , (22) 

(5^,)^;^'^ = --^£rfr'(xS)A^nd(x(r),x(r'))(C^).(r,T') . (23) 

Therefore the modification of energy shifts by the presence of the conducting plane boundary 
in a thermal bath is given by 

{SE,t:f = {SE.tf + {SE.tf . (24) 

III. MODIFICATION OF ENERGY SHIFTS BY THE PRESENCE OF A BOUND- 
ARY IN A THERMAL BATH 

In this section we will apply the formalism developed in the previous section to calculate 
separately the contributions of the thermal fluctuations and radiation reaction to the energy 
level shifts of a static two-level atom induced by the presence of a boundary in a thermal 
bath and the force acting on the atom as a result. We first give the basic results and then 
proceed to detailed discussions of the behaviors of the position-dependent energy shifts and 
thus force on the atom both in the low temperature limit (when the wavelength of the 
thermal photons is much larger than the transition wavelength of the atom) and the high 
temperature limit (when the wavelength of the thermal photons is much smaller than the 
transition wavelength of the atom). 
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A. Basic results 



For an atom at rest, we can compute the boundary-dependent two point function of the 
electric field in a thermal bath from the general form Eq. fl?I]) to get 

fc=— oo 

where u = t — r' . Consequently, the boundary-dependent statistical functions of the electric 
field can be expressed as 

he ^ / {5ij — 2ninj){uc + ikf^^ + Az'^ 



27r%o ^ V \(uc + ik/3-te)^-Az^P 

k=—oo 



{6ij - 2ninj){uc + ik[3) + Az 
^ [{uc + ikp + ief -Az'^Y 

he f {6ij — 2ninj)u'^c^ + Az^ {6ij — 2ninj)u^c^ + 



(26) 



27r%o V - ief - Az'^f [{uc + isY - Az"^]" 



(27) 



The linear susceptibility of the electromagnetic field, Eq. (!27ll . is independent of the tem- 
perature. Let us note that only the diagonal components {i = j = x, y, z) of the statistical 
functions are nonzero, and the components parallel to the conducting plane, i.e., the xx and 
yy components, are equal. Substituting the nonzero statistical functions of the atom and the 
electric field into the general formulas fl22|) and fl23l) . we can obtain the boundary-dependent 
contributions of thermal fluctuations and radiation reaction to the energy shifts of atomic 
levels. 



1^ /»oo 



fc=— oo j d 



. {I -2njnj){uc + ikl3f + Az^ {I - 2njnj){uc + ikl3f + Az^' 

[{uc + ik(5 -lef -Az^f ^ [{uc + ikf5 + iey -Az-^f ^ 



and 



{6E,)lf = ^^VK6|/i,(0)M)r / rfw(e™ + e-™^ 
Att eo ^ Jo 



(1 - 2n,nj)u''c^ + Az^ (1 - 2njnj)u\'' + 4z=^ ' 
^ [{uc - tey - Az^f [{uc + tey -Az^f ^ ' 



Here we have extend the range of integration to infinity for sufficiently long times t — tq. The 
integration can be evaluated using contour integrals with the help of the residue theorem. 
With a definition of the atomic static scalar polarizability 



the results can be cast into 



j j d 



2|(&|/^.(0)M)P 



(30) 



for the contribution of thermal fiuctuations to the energy level shift of the excited state 



(31) 



(5E_ 



, bnd 
'tf 



for that of the ground state, and 



1 + 



_ I ) /j(^o, z) - gji^o, z, (3) 



(32) 



ShujQaj 



(33) 



for the contribution of radiation reaction to the energy shifts of both the ground and excited 
states. Here we have defined 



fx{uJo,z) = fy{uJo,z) 



— cos(22:u;o/c) — sm(2zcJo/c) 

2 c z c 



(34) 



fz{i^o, z) = cos{2zLJo/c) - ^ sin(22:u;o/c) 



(35) 



gr,{uJo,z,p) = gy{uo 



[{uc + k(3y + Az^Y 



(36) 



gz{ujo,z,p) 



64c 



E 



du 



[{uc + k(5Y + Az'^Y 



(37) 



and summation over repeated indexes is implied. We note that the contribution of the 
radiation reaction to the energy shifts, Eq. fl33|) . does not depend on the temperature and 
the contribution of fluctuations of fields to the energy level shift of the ground state is just 
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opposite to that of the excited state. Adding up the contributions of thermal fluctuations 
and radiation reaction, we arrive at the position-dependent energy shift, 



1287r£n 



(38) 



for the excited state, and 



' 1287r£o 



(39) 



for the ground state. Since the energy shifts depend on the distance z between the atom and 
the conducting plane, the atom feels a force near the wall. In what follows, we examine in 
detail the behaviors of this force in various circumstances, considering separately the ground 
and excited atoms. 



B. Low temperature limit 

Let us now discuss the case when the wavelength of thermal photons is much larger than 
the transition wavelength of the atom, Aq, 

/5 » Ao , (40) 

where Aq = c/uq. This condition, which makes transitions from the ground state to the 
excited state virtually impossible, is very well satisfied at the usual room temperature. As 
a matter of fact, the first optical resonance of Rb atoms is as high as 1.8 x 10^ K. We will 
analyze how the interaction energy or the force acting on the atom behaves as the distance 
varies in three different regimes: the short distance, where the distance z is much smaller 
than the transition wavelength of the atom {z <S Aq), the intermediate distance, where the 
distance z is much larger than the transition wavelength of the atom but much smaller than 
the thermal wavelength (Aq ^-2^/9), and the long distance, where the distance z is much 
larger than the thermal wavelength {z ^ (3 ^ Aq). For these three distinct regimes, we will 
discuss the force on a ground state atom and an excited one separately. 
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1. Ground state atom 



Let us start with the short distance regime, i.e., when z ^ Aq. In this case, thermal 
fluctuations and radiation reaction contributions to the boundary- dependent energy shift 
become approximately. 



h 



AlT6r 



27r^c 



UJ, 



-^{a^ + ay- a^) log(2;u;o/c) 



+ 



15/54 63cu2/36 



[a^ \ ay - a^) + ———{2a^ + 2ay 



315/3 



6 



a. 



(41) 



and 



\bnd 
Jrr 



h ?>UJq 



(42) 



47r£o 32^3 

Here the contribution of radiation reaction is much larger than that of thermal fluctuations, 
and plays the dominant role in the total energy shift. Consider for example, a Rb atom 
{ujq = 2.37x lO^^s^^) placed 1 nanometer from the conducting boundary at room temperature 
(T = 300i^'), then the contribution of radiation reaction is about 200 times as large as that 
of thermal fluctuations. The closer the atom is to the boundary, the larger is this relative 
ratio. 

Since we are interested in the force on the atom, we will drop the 2— independent terms. 
Then we find by adding up the two contributions 



\bnd 
hot 



h 



Aixer 



?>UJq 

^(a. + a, + 2a.)-^^^^^ 

327r5cz2 ^ 
+ (2a^ + 2ay 



^(aa; + ay) log(2;u;o/c) 



315/5 



6 



(43) 



The first (leading) term on the righthand side of Eq. fH3l) is just the atom-wall interaction 

I 

energy at zero temperature in the short distance limit obtained in Ref. |2i] (an extra factor 
of h/ {4:7^60} is caused by the SI units used in the present paper) while the last term gives 
leading thermal corrections. The leading position dependent energy shift leads to the van 
der Waals force and the calculations above show clearly that the source of this force is purely 
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radiation reaction. The energy shift given in Eq. (H3!l induces an attractive force between 
a ground state atom and the plane wall. It is also interesting to note that the atom's 
polarization in the longitudinal and transverse directions are weighted differently in terms of 
their contribution to the force. For a ground state atom polarized isotropically the position 
dependent interaction energy shift reads 



{5E_f^f ^ ^^^^ ^ , g^^^g 



(44) 



where the first term is just the usual van der Waals result. One should bear in mind, however, 
that the polarization in different directions contributes differently in the van der Waals force 
even for isotropically polarized atoms. In fact, the van der Waals force on a ground state 
atom polarizable only in z direction would be twice as that polarizable only in the parallel 
directions. This is one of the new understandings we gain by not treating the atom as a 
limiting case of a macroscopic dielectric. 

Now let us turn to intermediate distance regime, i.e., when \q <^ z <^ j3. Here one finds 



(5E.)y 



h 



3uj^{a^ + ay) 3uJo 



32^3 



'a^ + + 2a^) ) cos(2zco'o/c 







+ fty + 2^2) sin(22;to'o/c) + 



3271^0^2 



3c 



and 



bnd 



h 



A-Ker 



315/3^ 



3ul{a^ + ay 



a. 



(45) 



3c^o 
32^3 



\oLx + ay + 2az) J cos(22;u;o/c) 



3uj 



\(a^ + ay + 2a^) sin(2zcJo/c) 



(46) 



Both the contributions of the thermal fluctuations and the radiation reaction have terms 
which are oscillating functions of z. However, these oscillating components cancel out when 
added up, leading to 



\bnd 
I tot 



h 



A-keq 

where the first term is comp 
action energy given in Ref. 



3c 



?>2-n^cz' 



(2a^. + 2ay - a^) 



(47) 



etely coincident with the zero temperature Casimir- Polder inter- 
2! and the atom's polarization in different direction contributes 
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equally. The second gives the thermal correction which is equal to that in the short dis- 
tance regime and it however depends on the atom's polarization direction. One sees that 
the Casimir-Polder force is a result of a residue energy shift due to vacuum fluctuations 
after cancelations of competing oscillating terms from thermal fluctuations and radiation 
reaction. This demonstrates that the Casimir-Polder force originates from a net result of 
vacuum fluctuations at zero temperature and radiation reaction. Here there is no clear dom- 
ination of one effect over the other as is in the short distance regime. In fact, the two effects 
are now comparable and opposing each other, and the net effect is much smaller than one 



either of them. We demonstrate this feature graphically in Fig. 1(a) and Fig. 1(b) For an 
isotropically polarized atom, we obtain 

h 



{SE_ 



\bnd 
hot 



3cao _|_ S2'K^caoz^ 



(4J 



which is always negative and the force on the ground state atom is again attractive in the 
intermediate distance regime. 

Finally, for the long distance regime, i.e., when z P ^ Xq, one has 



,bnd 
'tf 



h 



3u^{ax + ay 



ic^z 



[ax + ay + 2^2)^ cos{2zujo/c) 



16cz^ 



{ax + ay + 2az) sin(22;ci;o/c) + 



3c 



lQz^(3 



[ax + ay + 2az 



(49) 



and 



{6E. 



\bnd 



h 



3uj^{ax + ay) 



^(?Z 



3^0 
32^3 



{ax + + 2^2)^ cos(22;ci;o/c) 



3o;2 



IGcz'^ 



[a. 



+ ay + 2az) sin(22u;o/c) 



(50) 



Similar to the intermediate distance regime (Aq <C 2; <C /3), the oscillating terms in the 
contributions of thermal fluctuations and radiation reaction balance each other. So, again, 
in the present case, there is no clear domination of one effect over the other. However, the 
position-dependent energy shift of the ground state now has the form 

h 3c , „ s 1 3kBT 



A ir {ax + ay + 2a, 



{ax + ay + 2az) . 



(51) 
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This is linear in the temperature and depends on the polarization direction of the atom too. 
In the present case, the contributions from vacuum fluctuations are negligible as compared 
to that of thermal fluctuations in clear contrast to the intermediate regime. Therefore, in 
the low temperature limit {j3 ^ Aq), only when the distance between the ground state atom 
and the conducting boundary is much larger than the thermal wavelength {z ^ (3), do the 
thermal fluctuations become the dominating source of the position-dependent energy shifts 
and thus of the force acting on the atom. It is particularly interesting to note that the force 
acting on a ground state atom which is transversely polarized is twice as much as that of 
the atom which is longitudinally polarized. For an isotropically polarized ground state atom 
the total contributions 



This is in agreement with the Lifshitz's result [19|, which was obtained in a macroscopic 
approach to the problem where the force between two dielectric semi-infinite spaces at finite 
temperature was considered and the force on an atom is obtained as a limiting case when 
one dielectric is sufficiently rarefied. Note that the force is again attractive. 

2. Excited state atom 

After having analyzed the force on an atom in the ground state in the low temperature 
limit, let us now discuss what happens if the atom is in the excited state. One notes, from 
Eqs. (l3T])-(l33l). that the contributions of thermal fluctuations to the energy shifts for the 
ground and excited state atoms are of the same magnitude but opposite signs, while that of 
the radiation reaction are equal. Since radiation reaction is temperature independent, the 
thermal corrections to energy shift of the excited state come solely from thermal fluctuations 
and thus will be just opposite to that of the ground state. As a result, in the short distance 
regime where the contributions from radiation reaction dominate over thermal fluctuations 
(and vacuum fluctuations), one finds 



22^3 [ax + ay + 2a^) 3x5^6 ^^'^^ + '^'^v ~ 

15 



(53) 



which agrees with the result for the ground state atom in the leading order (refer to Eq. (H3|) ). 
Therefore, in the low temperature and short distance limit, the force acting on the atom 
reduces to the Van der Waals result for both the ground and excited in the leading order 
plus a small thermal corrections which differs in sign. The net force is attractive. It is 
worth pointing out here that the force on excited atoms, in general, may not be necessarily 
the gradient of a potential. However, in the intermediate distance regime, where radiation 
reaction no longer dominates, the total energy level shift deviates considerably from that of 
the ground state as a result of the fact that the oscillating terms from thermal fluctuations 
and radiation reaction do not cancel as they do in the case of the ground state atom. The 
result is 



bnd 
tot 



h 



'd>ujl{ax + ay) 3t^o 
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Amplitude of the oscillating terms is much larger than the Casimir-Polder term which be- 
haves like Xj and consequently, for a given atom, the position-dependent energy shift of the 
excited state may be much larger than that of the ground state depending on the distance z 
of the atom from the boundary. The thermal corrections are very small compared with the 
vacuum fluctuation contributions. For an isotropically polarized atom, we can write Eq. (!54l) 
as 

h 



bnd 
tot 



Scan 327r^cao2^ 
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(55) 



which may be positive or negative depending on the distance z so that the force acting on the 
excited state atom may be attractive or repulsive and even be zero. This is quite different 
from the case of a ground state atom, where the position-dependent energy shift is always 
negative and the force is attractive. It is interesting to note that, for a given atom, there 
exist certain values of the distance z such that the oscillatory terms sum up to zero and 
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when that happens the position-dependent energy shift becomes 
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which leads to a force on the excited atom equal in magnitude to that on ground state atom 
but opposite in direction. It is remarkable that whenever the oscillatory terms do not sum 
up to zero, the force on the excited atom is generally much stronger than the usual Casimir- 
Polder force on the ground state atom and alternates between attractive and repulsive as 
the distance varies. The total energy shift for an excited state is plotted in Fig. 
with Fig. [1] to show this. 
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Now it is time to look at the long distance regime. It is easy to show that the energy shift 
can be approximated as 
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Here the oscillating terms arising from vacuum fluctuations are the same as in the interme- 
diate distance regime but the thermal correction term {l/z^f3) overtakes the Casimir- Polder 
one {l/z"^) . A comparison with Eq. (IST]) shows that, although the thermal effects plays a 
major role in the position-dependent energy shift for the ground state, it is generally not the 
case for the excited state because of the presence of oscillating terms. In fact, the thermal 
correction term is of higher order when compared with the oscillating terms, unless the atom 
is such placed that the oscillatory terms sum up to zero. Therefore, the net contributions 
of vacuum fluctuations and radiation reaction in general dominate even in the long distance 
regime if the atom is in the excited state, as opposed to the ground state where thermal 
effects dominate in the long distance regime. For an isotropically polarized atom. 
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which again can be positive or negative depending on the position of the atom. 

The above discussions reveal that the force acting on the excited atom is definitely at- 
tractive only in the short distance regime and can be either attractive or repulsive in the 
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intermediate and long distance regime. This is in a clear contrast to the ground state atom 
where the force is attractive in all regimes. When the oscillatory terms sum up to zero, we 
have 



\bnd h cao 1 aoksT 



leading to a force which is equal to the Lifshitz result in magnitude but opposite in direction. 



C. High temperature limit 

Now we turn our attention to the high temperature limit, where we assume the wavelength 
of the thermal photons to be much smaller than the transition wavelength of the atom, i.e., 

< Ao , (60) 

which can also be written as I3ujq/c <^ 1. Then Eqs. fl3Tl) and fl32|) . for the contributions of 
the thermal fluctuations to the position-dependent energy shifts, can be approximated as 

2c 
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As in the low temperature limit, we will separately examine the position-dependent energy 
shifts and thus the force acting on the atom for both the ground and excited states in three 
different distance regimes. 

We first consider the case when the atom is placed so close to the boundary that the 
distance z is much smaller than the thermal wavelength {z <^ (3 <^ Aq). It follows that 
contributions of thermal fluctuations and radiation reaction to the energy level shifts are 
now respectively. 



and 



{6E^tf = {SE^tf ^ + ^y + 2«^) • (63) 
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Here the contributions of thermal fluctuations (l62l) are much larger than that of the radiation 
reaction fl63l) and play the dominating role in the position-dependent energy shifts for both 
the ground and excited state. It is interesting to note that in the high temperature limit, even 
when the distance z is much smaller than the wavelength of thermal photons, the thermal 
effects dominate over zero-point contributions. This is in contrast to the low temperature 
limit where thermal effects dominate only in the long distance regime when the distance z 
is much larger than the wavelength of thermal photons. For an isotropically polarized atom, 
we obtain 

= -i5E^\ot - 4^4^ = ^^^^ ' (64) 

which is linear in temperature, and gives a positive position-dependent energy shift for a 
ground state atom and a negative one for an excited state. This means that the force acting 
on an atom is repulsive if the atom is in the ground state and attractive if otherwise. Notice 
that the force on a ground state atom in the present case is different from the Lifshitz's 
result in the low temperature and long distance limit in direction although it is equal in 
magnitude. 

If the distance z is much larger than the thermal wavelength but much smaller than the 
transition wavelength, i.e., /5 <C z <^ Aq, then we find 
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For the total position-dependent energy shifts, we can obtain 
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for the ground state, and 
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for the excited state. Here there is again no strong domination of one effect over the other. 
However, the above resuhs show that now whether the force acting on the atom is attractive 
or repulsive depends crucially on the polarizations of the atom. For an example, a change 
of relative ratio of polarizations in different directions may result in a change of direction of 
the force. We note that for an isotropically polarized atom, the first term in both Eg. (1671) 
and (168|) is zero, so 
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Here for an excited state atom, the position-dependent energy shift is always negative and 
induces an attractive force, but for a ground state atom, it is determined by the magnitude 
of the distance z and the wavelength of the thermal photon. The crucial dependence of the 
direction of the force acting the atom on its polarization in this case is an interesting feature 
of atom-wall interactions revealed in the present paper. 

Thirdly, for the long distance limit such that the distance z is much bigger than the 
transition wavelength {z ^ Xq ^ (3), the contributions of thermal fluctuations and radiation 
reaction can be written as 
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Clearly, in the present case, the contribution of radiation reaction is much smaller than that 
of thermal fluctuations since [3ujq/c <^ 1 and thus the vacuum fluctuation effect dominates. 
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So the position-dependent energy shifts of both the ground and excited states are dominated 
by the contribution of thermal fluctuations. Assuming isotropic polarization, we arrive at 
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(73) 



which has both terms of oscillating function of distance z and a Lifshitz-like term. If the 
atom is so placed that the oscillating part is equal to zero, the position-dependent energy 
shift of the ground state reduces to the Lifshitz's result. 

In all the discussions above, we have only studied energy level shifts of an atom in its 
ground and excited state separately. An interesting issue yet to be addressed is the thermal 
average of energy level shifts of an atom in the equilibrium with thermal photons. This is 
particularly desirable in the high temperature limit, since now the ground state atoms have 
a high possibility to absorb thermal photons to transit to the excited state. Taking into 
account the fermionic nature of the two-level atom, which can be seen by noting that the 
atomic raising and lowering operators obey the anticommutation relation, {5+, S-} = 1, we 
can calculate the thermal average of contributions of the thermal fluctuations and radiation 
reaction to the position-dependent energy shifts as follows. 
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When the temperature is low, the thermal average shift is essentially the shift of the ground 
state as expected, and significant deviations only occur when temperature is high. So, in 
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what follows, we will examine in detail the average energy shifts in the three distinct distance 
regimes in the high temperature limit. 

First, the short distance regime, where the distance z is assumed to be much smaller than 
the thermal wavelength, i.e., -C /5 <^ Aq. Now the contributions of thermal fluctuations 
and radiation reaction to the average energy shifts can be approximated as 

^ J_ _ (a, + a, + 2a,) , (76) 

^ '^^ 47r£oV32^3 64czV ^ ^ ' ^ ' 

and 

Here one can see that the contribution of temperature-independent vacuum fluctuations to 
the average energy shift cancels that of the radiation reaction when added up, leaving a 
correction due to thermal fluctuations as the dominant term in the average energy shift. So 

This means that the average energy shift decreases as the temperature increases in the short 
distance regime (z -C /3 <^ Aq). For isotropic polarization, Eq. ( 1781) becomes 

This is always negative and leads to an attractive average force on the atom. 

Second, we consider the intermediate distance regime, where the distance z is much larger 
than the thermal wavelength but much smaller than the transition wavelength of the atom, 
i.e., /3 -C 2; -C Aq. Now one finds 
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so the average energy shift is given by 
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Here the average energy shift is caused mainly by the radiation reaction and the thermal 
correction due to thermal fluctuations is quadratic in (3 in contrast to the linear dependence 
in the short distance regime. For an isotropically polarized atom, Eq. (15^ reduces to 



(83) 



where the first term is just the Van der Waals result. 

Finally, we can estimate the contributions of thermal fluctuations and radiation reaction 
to the average energy shift, in the long distance regime (2; ^ Aq ^ /?), to get 
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Noting that the oscillating terms in the contributions of thermal fluctuations are canceled 
out by that of the radiation reaction in the calculation of the total shift, we obtain 
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This coincides with Eq. (!82|) for the case of /? ^ z <^ Aq in mathematical appearance but 
differs in physical origin. The reason is that here the average energy shift comes mainly from 
the contribution of the thermal fluctuations while in the intermediate distance regime the 
radiation reaction plays the major role. 



IV. CONCLUSION 



Using the DDC formalism which separates the contributions of thermal fluctuations and 
radiation reaction and allows a distinct microscopic treatment to atoms in the ground and 
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excited states, in contrast to the macroscopic approach where atoms are treated as a hmiting 
case of a dielectric, we have investigated the modification by the presence of a plane wall 
of energy level shifts of two-level atoms which are in multipolar coupling with quantized 
electromagnetic fields at the finite temperature. The position dependent energy shifts give 
rise to an induced force acting on the atoms. We have examined the behaviors of this force 
in both the low temperature limit and the high temperature limit for both the ground state 
and excited state atoms. 

In the low temperature limit where the wavelength of thermal photons is assumed to be 
much larger than the transition wavelength of the atom, our calculations reproduce, for a 
ground state atom, the van der Waals force at short distance regime, the Casmir-Polder 
force in the intermediate regime and the Lifshitz result in the long distance regime, and 
show that the force is attractive in all three regimes. However if the atom is in the excited 
state, the atom-wall force agrees with that in the ground state in the leading order only 
in the short distance limit. In the intermediate distance regime the force on the excited 
state atoms is generally strong than the Casimir-Polder force on the ground state atoms and 
alternates between attractive and repulsive as the atom's distance from the wall varies. In 
the long distance regime, the main source of the force is the net contributions of competing 
vacuum fluctuation and radiation reaction and the thermal correction term is in general of 
higher order. This is in clear contrast to the ground state atom where the thermal effect 
(the Lifshitz term) dominates. Furthermore, the force can again be attractive and repulsive 
depending on the atom's distance. 

In the high temperature limit where the wavelength of thermal photons is assumed to be 
much smaller than the transition wavelength of the atom, we find that the thermal effects 
dominate over zero-point contributions even in the short distance limit where the distance is 
small as compared to the transition wavelength and the force is linear in temperature in the 
leading order in all three regimes. In the short distance regime the force on an excited atom 
is attractive while that on a ground state atom is repulsive. In the intermediate distance, 
the direction of the force acting on the atoms depends crucially the polarization. Only in 
the long distance regime does the force have oscillatory part, leading to a force which would 
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either be attractive or repulsive. 

We also have calculated the thermal average of energy level shifts for an ensemble of 
atoms in equilibrium with the thermal photons in the high temperature limit. We find 
that the average force is attractive and decreases as 1/T as temperature T grows in the 
short distance regime. However, both in the intermediate and long distance regimes, the 
average force coincides with the van de Waals force in the leading order, which is curiously 
independent of the temperature. One should note, however, that although the force has the 
same form in two different regimes, the physical origin is quite different, with the force in 
the intermediate regime being generated by radiation reaction and that in long distance by 
thermal fluctuations. 

By not treating atoms as a limiting case of a macroscopic dielectric, we are able to reveal 
the influence of different atomic polarizations on the energy level shifts and thus the force 
acting on the atoms. Interestingly, we find that both the magnitude and the direction of the 
force acting on an atom have a clear dependence on the atomic polarization directions. In 
certain cases, a change of relative ratio of polarizations in different directions may result in 
a change of direction of the force. 

Finally, as far as the separation of effects of fluctuations and radiation reaction in the 
DDC prescription itself is concerned, we only flnd a clear domination of one over the other 
in short distance regime in the low temperature limit where radiation reaction effect domi- 
nates and in both the short and long distance regimes in the high temperature limit where 
fluctuations effect does. An interesting issue for future discussions is the non-equilibrium 
situation where the atom and fleld have different temperatures and the two terms can thus 
be varied independently. In fact, in this case, the contribution of radiation reaction (atom 
fluctuations) would only depend on the atomic temperature and the thermal fluctuation 
contribution (fleld fluctuations) would only depend on the fleld temperature. We hope to 
return to this issue in the future. 
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FIG. 1: Energy level shifts of both a ground state atom and an excited one polarized isotropically in 
the intermediate distance regime (Aq <C 2; <C /?) at room temperature T = 300if. Here the frequency 
of Rb atom, loq = 2.37 x lO^'^^s^"'^ is used, and the levels shifts are in the units of ao/(4vreo)- (a) 
The contributions of thermal fluctuations (dash-dotted line) and radiation reaction (dashed line) 
to the level shift of a ground state atom is plotted, demonstrating the magnitude of two effects are 
comparable and there is no clear domination of one over the other, (b) The level shift of a ground 
state atom as a net result of two competing effects is plotted. The net result is clearly much less 
than any of the competing effects, (c) The energy level shift of an excited state atom, Eq. (jSSp . is 
plotted, showing clearly differences both in magnitude and sign from that of a ground state. 
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